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Abstract. We construct a continuous semigroup of weak, dissipative solutions to a nonlinear 
partial differential equations modeling nematic liquid crystals. A new distance functional, deter- 
mined by a problem of optimal transportation, yields sharp estimates on the continuity of solutions 
with respect to the initial data. 
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In this paper we investigate the Cauchy problem 



oo J X 



Formally differentiating the above equation with respect to the spatial variable x, we obtain 

s 



(lit + uu x ) x = - ul, (1.2) 



X 

$_i ' whereas yet another differentiation leads to 



^txx 2'U'x'U'xx H'H'xxx — 0. (1 .3) 

Either of the forms (1.2) and (1.3) of the equation in (1.1) is known as the Hunter-Saxton equation. 
In this paper we analyse various concepts of solutions for the above equations, and construct a 
semigroup of globally defined solutions. Moreover, we introduce a new distance functional, related 
to a problem of optimal transportation, which monitors the continuous dependence of solutions on 
the initial data. 

Physical significance 

The Hunter-Saxton equation describes the propagation of waves in a massive director field of 
a nematic liquid cristal [HS], with the orientation of the molecules described by the field of unit 



1 



vectors n(t,x) = (cos u(t,x), sin u(t,x)^j, x being the space variable in a reference frame moving 
with the linearized wave velocity, and t being a slow time variable. The liquid crystal state is a 
distinct phase of matter observed between the solid and liquid states. More specifically, liquids 
are isotropic (that is, with no directional order) and without a positional order of their molecules, 
whereas the molecules in solids are constrained to point only in certain directions and to be only in 
certain positions with respect to each other. The liquid crystal phase exists between the solid and 
the liquid phase - the molecules in a liquid crystal do not exhibit any positional order, but they do 
possess a certain degree of orientational order. Not all substances can have a liquid crystal phase 
e.g. water molecules melt directly from solid crystalline ice to liquid water. Liquid crystals are 
fluids made up of long rigid molecules, with an average orientation that specifies the local direction 
of the medium. Their orientation is described macroscopically by a field of unit vectors n(i, x)). 
There are many types of liquid crystals, depending upon the amount of order in the material. 
Nematic liquid crystals are invariant under the transformation n i— > — n, in which case n is called a 
director field, so that the rodlike molecules have no positional order but tend to point in the same 
direction (along the director). The director field does not remain the same but generally fluctuates. 
Obtaining the equation governing the director field represents the crucial point for the modeling 
of nematic liquid crystals since it is advantageous to study the dynamics of director field instead 
of studying the dynamics of all the molecules. The fluctuations of the director field are mainly 
due to the thermodynamical force caused by elastic deformations in the form of twisting, bending, 
and splay (the latter being a fan-shaped spreading out from the original direction, bending being a 
change of direction, while twisting corresponds to a rotation of the direction in planes orthogonal 
to the axis of rotation). Consider director fields that lie on a circle and depend on a single 
spatial variable x so that twisting is not allowed. To describe the dynamics of the director field 
independently of the coupling with the fluid flow, let u(t, x) be the perturbation about a constant 
value. The asymptotic equation for weakly nonlinear unidirectional waves is precisely equation 
(1.2), obtained as the Euler-Lagrange equation of the variational principle 

5 / (u t u x + uu 2 x ) dxdt = 0, 
Jti J R 

for the internal stored energy of deformation of the director field if dissipative effects are neglected 
(corresponding to the case when inertia effects dominate viscosity) - see [HS] for the details of the 
derivation. Unlike other studies, in the Hunter-Saxton model the kinetic energy of the director 
field is not neglected. In the asymptotic regime in which (1.2) is derived, the nondimensionalized 
kinetic energy density is u x so that the condition 

u x (t, x) dx < oo (1-4) 

has to hold at any fixed time t for a physically meaningful solution to the Hunter-Saxton equation. 

Equation (1.1) is also relevant in other physical situations, e.g. it is a high-frequency limit of 
the Camassa-Holm equation [DP], a nonlinear shallow water equation [CH, J] modeling solitons 
[CH] as well as breaking waves [CE] . 

Geometric interpretation 

An interesting aspect of the Hunter-Saxton equation (see [KM]) is the fact that, for spatially 
periodic functions, it describes geodesic flow on the homogeneous space Diff(<S)/Rot(<S) of the 
infinite-dimensional Lie group Diff(<S) of smooth orientation-preserving diffeomorphisms of the 
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unit circle S modulo the rotations Rot(<S), with respect to the right-invariant homogeneous metric 

(/> 9) = J fxQx dx. The geometric interpretation of the Hunter-Saxton equation establishes a 

natural connection with the Camassa-Holm equation, which describes geodesic flow on Diff(<S) with 

respect to the right-invariant metric (/, g) = J (fg + f x g x ) dx, see [K, CK]. A similar geometric 

interpretation of (1.1) on the diffeomorphism group of the line holds also for smooth initial data u 
in certain weighted function spaces but the involved technicalities are more intricated (see [C] for 
the case of the Camassa-Holm equation). 

Integrable structure 

The Hunter-Saxton equation has a an integrable structure. The equation has a reduction (see 
[BSS, HZ1]) to a finite dimensional completely integrable Hamiltonian system whose phase space 
consists of piecewise linear solutions of the form 

n 

(t,x) = ^2ai(t)\x-Xi(t)\, (1.5) 



with the constraint 



the Hamiltonian being 



J>i(t) = 0, (1.6) 



i=i 



1 n 

H(x,a) = - £ 
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Due to their lack of regularity, functions of the form (1.5) are not classical solutions of (1.2). 
Below we will discuss in what sense they are weak solutions of the Hunter-Saxton equation. Let 
us point out that the constraint (1.6) is the necessary and sufficient condition to ensure that the 
distributional derivative x i-> u x (t,x) of a function of the form (1.5) belongs to the space L 2 (M). 
Thus (1.4) holds. 

In the family of smooth functions u : M t— > M all of whose derivatives d™u decay rapidly 
as x — > ±oo, the Hunter-Saxton equation is bi-Hamiltonian [HZ1]. If D~ x is the skew-adjoint 
anti-derivative operator given by 

(D- 1 /) (x) = \[f_ ~ ^°°) f(x) dx, f E V(R), 

the first Hamiltonian form for the Hunter-Saxton equation is 

S1~L 1 f 

u t = J\ - r ^-, Ji = u x D~ 2 - D~ 2 u x , H\ = - I u 2 x dx, 
du 2 J R 

whereas the second, compatible, Hamiltonian structure is 



Ut = J2^^-, J 2 = D 1 , 7i 2 = ^ I uu 2 



X 



Moreover, the Hunter-Saxton equation is formally integrable e.g. it has an associated Lax pair 
(see [BSS]). However, the complete integrability of the equation has been established only in the 
previously mentioned case when it reduces to a finite dimensional dynamical system. 
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The notion of solution 



Physically relevant solutions of the Hunter-Saxton equation need to be of finite kinetic energy 
so that (1.4) must hold. This leads naturally to functions u(t,x) with distributional derivative 
u x (t, •) square integrable at every instant t. Note that the integrability assumption u x (t, •) G L 2 (M) 
already imposes a certain degree of regularity on the function u. This suggests that it might be 
possible to incorporate a reasonably high degree of regularity in the concept of weak solution to the 
Hunter-Saxton equation. Let us first consider the concept of weak solutions introduced by Hunter 
and Zheng [HZ2] . 

Definition 1.1 A function u(t,x) defined on [0, T] x iR is a solution of the equation (1.2) if 

(i) u G C([0,T] x M;M) and u(0, x) = u(x) pointwise on M; 

(ii) For each t G [0, T], the map x u(t,x) is absolutely continuous with u x (t,-) 6 L 2 (1R). 
Moreover, the map t u x (t,-) belongs to the space L°°([0, T]; L 2 (M)) and is locally Lipschitz 
continuous on [0, T] with values in H^ C (M); 

(Hi) Equation (1.2) holds in the sense of distributions. 

Here and below, by a mapping / that is locally Lipschitz or locally bounded on [0, T] with values 
in Hj^ c (iR) we understand the following: for every n > 1 there is a constant K n > such that 



sup 

WeT>(-n,n): || </>||„i (JR) <1} 



(f(t)-f(s),i>) <K n \t-s\, t,se[0,T], 



respectively 



sup 

WeV(-n,n): || (M) <1} 



</(*), V) <K n , t€[0,T\. 



Here V(a, b) is the family of smooth functions / : 1R — ► 1R with compact support within (a, b) C M. 

To a function u : [0,T] x 1R — > iR with the above properties associate the function F : 
[0, T] x iR -» M defined by 

F & x ) = \(f ~ J™) U * dx - (L7) 

Then F G L^ c ([0,T] x M;M) C Lf oc ([0,T] x iR;iR). Moreover, = | u\ so that equation (1.2) 
becomes 

(u t + uu x - F) x = (1.8) 

in the sense of distributions. Note that uu x G L 2 oc ([0, T] x iR; iR). From (1.8) we infer the existence 
of a distribution h(t) so that u t +uu x — F = h(t)®l(x), where l(x) stands for the constant function 
with value 1 on iR. If H(t) is a primitive of the distribution h(t), we deduce that the distribution 
U = u - H(t) <g) l(x) satisfies U t = u t - h(t) (g) l(x) = F - uu x G lf oc ([0,T] x iR; iR) and 
U x = u x e L 2 oc ([0,T] x R;1R). Therefore C7 G if/ oc ([0,T] x iR). Moreover, since U t = F - uu x G 
^7oc([0' -Tl' ^loc (-^)) ensures that [/ is locally Lipschitz as a function from [0, T] to H^ o l c {lR) and so 
is also n, we deduce that h(t) l(x) = u — U shares this property too. But then h : [0, T\ — > iR has 
to be Lipschitz continuous. We infer that u = U + H(t)(g>l(x) belongs to the space H} oc {[Q, T] x iR). 
Since the requirement (iii) in Definition 1.1 ensures that the identity 



J J (4>xtu + i 4> xx u 2 - i (fmfj dxdt = 
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holds for every smooth function cf> : (0, T) x M — > M with compact support in (0, T) x M, we see 
that the notion of weak solution in the sense of Definition 1.1 is stronger than the concept of weak 
solution introduced by Hunter and Saxton [HS]. Another useful conclusion that can be drawn from 
the previous considerations is that for a function u with regularity properties specified in (i)-(ii) of 
Definition 1.1, the requirement (hi) from Definition 1.1 is equivalent to asking that the equation 

u t + uu x = F + h(t) o l(x) (1.9) 

holds in distribution sense for some Lipschitz continuous function h : [0,T] — > M. Any such 
function h is admissible. Among all these possibilities the most natural one corresponds to the 
special choice h = 0. This leads us to the form (1.1) of the Hunter-Saxton equation. 

In the following, we say that a map 1 1— > u(t, •) from [0, T] into Lf oc (M) is absolutely continuous 
if, for every bounded interval [a, b], the restriction of u to [a, b] is absolutely continuous as a map 
with values in L p ([a, b]) . We can thus adopt the following notion of a weak solution. 

Definition 1.2 A function u(t,x) defined on [0, T] x M is a solution of the equation (1.2) if 

(i) u 6 C([0, T] x M;M) and u(0,x) = u(x) pointwise on M; 

(ii) For each t £ [0, T], the map x i— > u(t,x) is absolutely continuous with u x (t,-) 6 L 2 (1R). 
Moreover, the map 1 1— > u x (t, •) belongs to the space L°°([0, T]; L 2 (M)); 

(Hi) The map t u(t,-) 6 Lf oc (M) is absolutely continuous and satisfies the equation (1.1) 
for a.e. t 6 [0,T]. 

The concept of solution introduced in Definition 1.2 is stronger than that corresponding to Defini- 
tion 1.1. Indeed, for a function u satisfying all the requirements of Definition 1.2 we infer by (1.1) 
that u tx £ L^ c ([0,T};Hf o l(]R)) since u t = -uu x + F and uu x , F £ Lf£ c ( [0, T] ; Lf oc (M) ) . This 
yields that the map t i— > u x (t, •) is locally Lipschitz continuous on [0,T] with values in H^l(M). 
We thus recover the apparently missing part from the requirement (ii) in Definition 1.1. 

We remark that, even with this stronger definition, solution are far from unique. For example, 
consider the initial data 

u(x) = 0. (1.10) 
There are now two ways to prolong the solution for times t > 0. On one hand, we can define 

u(t,x) = x£M, t>0. (1.11) 

On the other hand, the function 

( -2t if x < -t 2 



u(t, x) 



2x_ 

t 



if Id < t 2 



for t > 



(1.12) 



2t 



if x>t 2 



provides yet another solution. To distinguish between these two solutions, we need to consider the 
evolution equation satisfied by the "energy density" u 2 ,, namely 



(u 2 x ) t + (uu 2 x ) x = 0. 



(1.13) 



For smooth solution, the conservation law (1.13) is satisfied pointwise. Notice that the solution 
defined by (1.10), (1.12) satisfies the additional conservation law (1.13) in distributional sense, i.e. 



/ / { u x <£t + uu x (fx] dxdt = 

J Jm+xR 



(1.14) 
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for every test function <p G C\ (iR+ x JR) whose compact support is contained in the half plane where 
t > 0. On the contrary, the solution defined by (l.lO)-(l.ll) dissipates energy. More precisely, for 
every t 2 > ti > we have 

/ u 2 (t 2 , x) (f(t 2 , x) dx — \ u x (t 2 , x) <p(t 2 , x) dx < / / {u 2 x p> t + uu 2 cp x } dxdt , (1-15) 

JR JR Jt! JR 

for every test function ip G Cl(M+ x M). In the sequel, we say that a solution is dissipative if the 
inequality (1.15) holds for every t 2 > t\ > 0, ip 6 Cl(lR + x 1R). Notice that the solution (1.10), 
(1.12) does not satisfy (1.15) when ti = 2, t 2 > 0. 

At this point in the discussion it is worthwhile to point out that the most important feature in 
the definition of weak solutions is the requirement (1.4). A continuous function u : [0, T] x M — > M 
with square integrable distributional derivative u x (t,-) belonging to the space L°°([0, T]; L 2 (]R)) 
is not necessarily bounded, nor does it have a pre-determined asymptotic behaviour at infinity, as 
one can see from the example 



u(t, x) 



\x\ 5 sin(|x| s ) if t > 0, \x\ > 1 , 

\x\i sin(\x\i) if t > 0, \x\ < 1 . 



Nevertheless, the possibility that some additional structural information about the behaviour of 
such functions at infinity might be inferred from some invariance properties of the Hunter-Saxton 
equation should be ruled out. To do this, consider solutions of the type (1.5) with the constraint 
(1.6). This type of solutions enter into the framework of Definition 1.2 and for any N(t) > 
max{|xi(£)|, |x n (£)|} we have 

u t (t,x) = F(x) a.e. on \x\ > N(t), 

so that for all j > N(t), 

Ut (t,j)-u t (-j) = F(j)-F(-j) = X - J[u 2 x (t,x)dx = X - j^u 2 x (t,x)dx. (1.16) 

But the quantity I = - f u 2 (t, x) dx is an invariant (time-independent) cf. [HZ1, BSS]. Moreover, 
2 Jr 

the special form of the solutions guarantees that at every fixed t > 0, 

n 

Uoo(t) = nm u(t,x) = y a.i(t)xi{t) = — lim u(t,x). 

x—>oo — J x—> — oo 

i=l 

and u^t) = u(t,j) = -u(t, -j) for all j > N(t). Thus (1.16) yields 

«oo(*) =«oo(0) + i/t, t>0. (1.17) 

Unless / = 0, in which case u is constant, we see from (1.17) that the asymptotic behaviour of the 
solutions changes with time. On the basis of this set of examples we conclude that the asymptotic 
behaviour of the solutions at infinity should not be prescribed a priori. However, the previous 
set of examples indicates that a possible restriction would be to require u 6 L°°([0, T] x M) if 
u 6 L°°(iR). In this case the space of functions introduced in Definition 1.2 (that is, bounded 
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functions with all the properties specified in Definition 1.2 except the condition that u satisfies 
equation (1.1) in L 2 [—n,n] for every n > 1) is a Banach space when endowed with the norm 



\U T 



sup {\u(t, x)\} + ess-sup / u x (t,x)dx. 
(t,x)e[o,T]xR te[o,T] Jr 



(1.18) 



It is also worth noticing that a function entering the framework of Definition 1.2 has further 
reqularity properties that are not explicitely stated. For example, we have the Holder continuity 
property 

\u(t, x) - u(t, y)\ < K(t) y/\x-y\, t>0, x,y € M, 
with K(t) = \\u x (t, -)||l2(jr), since 

\u(t,x)-u(t,y)\ 2 = [ V u x (t,()d( 2 < \x-y\ • F u\ (t, C) d( < \x - y\ [ u 2 x (tX)d(. 



2 - Global existence of dissipative solutions 

For twice continuously differentiable solutions, the derivative v = u x of the solution u to (1.1) 
satisfies the equations 

v 2 

v t + uv x = — -, (2.1) 

(v 2 ) t + (uv 2 ) x = 0. (2.2) 
Define the characteristic 1 1— > £(t,y) as the solution to the O.D.E. 
d 

- e(t, y) = u(t, at, y)) , e(0, y) = y. (2.3) 

From (1.2) it follows that the evolution of the gradient u x along each characteristic is described by 

± Ux (t,Z(t,y)) =-\u 2 x (t,at,y)). (2.4) 
Observe that the solution of the O.D.E. 

z = -z 2 /2, z{0) = z 

is given by 

- £k (2 - 5) 

If Zq > this solution is defined for alH > 0, whereas if zq < 0, this solution approaches — oo at 
the blow-up time 

T(z ) = -2/zo (2.6) 

Note that if u(x) ^ then there is some x G M with u(xo) < so that the characteristic cruve 
t l— ► £(t,u(xo)) will blowup in finite time. Nevertheless, if liminf ;re jR{^(a ; )} > — oo, then T > 0, 
where 

T = inf \^-}>0, (2.7) 

{xeiR. u x (x)<o} lu x (x)i 
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and on the time interval [0, To) the method of characteristics can be used to construct the unique 
solution of (1.1). Let us describe the construction in detail. From (2.3) we get 



£«- = «■"■« = (2 ' 8) 

since 

in view of (2.4) and the solution formula (2.5). The unique solution of the linear O.D.E. (2.7) with 
initial data £, x (0,y) = 1 is given by 

Ut,y)=(} + \u x (y)) 2 . (2.10) 

Since 1 + \ u x (y) > for t € [0, To), relation (2.10) shows that for each t 6 [0, To) the map 
y I— > £(t, y) is an absolutely continuous increasing diffeomorphism of the line. Define the absolutely 
continuous function ip by 

<p(v) = -j I sign(y-x)u 2 x (x)dx (2.11) 

4 J Ft 

so that <p x (y) = \u 2 x (y). Note that by (2.10), 

u 2 

tt x = u x + t^. (2.12) 
Since £ t (0, y) = as £(0, y) = y, integration of (2.12) with respect to the spatial variable x yields 

Zt(t,y) =u{y)+ t - I sign(y - x) u 2 x {x) dx (2.13) 

4 J Ft 

and thus 

f* t 2 f 

€(t,y) = y+ £ t (s,y)ds = y + tu(y) + — / sign(y - x) u x (x) dx. (2.14) 

JO 4 J Ft 

The value of the solution u along the characteristic curve 1 1— ► £(£, y) is 

u(t,£(t,y)) =u(y) + I / sign(y-x)u 2 x (x)dx. (2.15) 

1 J Ft 

This relation is obtained by combining (2.14) with (2.3). The increasing diffeomorphism of the 
line y i-> £(t,y) given by (2.14) and formula (2.15) yield the unique solution of the Hunter-Saxton 
equation on the time interval [0, To). The above approach works as long as 2 + t u x (x) > but 
breaks down at T = T with T given by (2.7). The reason for the breakdown is that 

liminfjinf u x (t,x)} = -co (2-16) 

£|To x€:Ft 

in view of (2.9) and the definition (2.7) of Tq. Note that at t = Tq we have might have £ x (t, x) = 
for all x € (a,b) C M so that the map y ^— > £(i, y) is not any more an increasing diffeomorphism of 
the line. Nevertheless, the previous considerations suggest the following approach in the general 
case when u x 6 L 2 (M), covering situations when possibly T = as it is the case for e.g. u(x) = 
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X3 (1 — x) 3 X[o,i] • Here xa stands for the characteristic function of the set A, defined by xa{x) = 1 
if x G A and xa{x) = if x A. 

Let u G C(M) be such that its distributional derivative u x is square integrable. Define <p : 
iR + x iR -> iR by 

Py(*»i/) = 2 S x(j/)X[u B >-2/t](l/) (2-17) 

so that 

<p(t,y) = - sign(y-x)u 2 x (x)dx, t > 0, (2.18) 

4 ./[u x >-2/t] 

with the understanding that 



¥>(0> y) = t / si g n (y ~ dx. 

4 JJR 



In other words, if u(t,£(t,xo) blows up before to > 0, then the point xq is not included in the 
domain of the integral defining tp(to,-), because 

T(u x (x)) > t if and only if u x (x) > —2/t , 

according to (2.4) and (2.6). Observe that (2.18) and Young's inequality yield 

\\<p(t,-)h-(R) < 7 / ul(x)dx, t>0. (2.19) 
4 Jm 

In the (t, x)-plane, the characteristic curve starting at y is obtained as 

£(t, y)=y + tu(y) + [ (t-s) <p(s, y) ds . (2.20) 
Jo 

The value of the solution u along this curve is 



u( 



(t, £(t,y)) = u{y) + [ <p{s,y)ds. (2.21) 
J o 

Observe that for alH > and y G iR, 

Ht(t,y) = u(t,at,y)) (2.22) 

in view of (2.20)-(2.21). 

Theorem 1. Given any absolutely continuous function u : M — > iR with derivative u x G 
L 2 (M), the formulas (2. 18)- (2. 20) provide a dissipative solution to (1.1), defined for all times 
t > 0. 

Proof. We proceed in several steps. First of all, for any fixed t > 0, the map y t— > £(t,y) is 
absolutely continuous since ip y (t, •) G L 2 (M). We claim that for any fixed t > the map y i-> £(£, y) 
is nondecreasing on iR with lim^-i-oo y) = ±oo. 

Indeed, if u x (y) > — | then u x (y) > — | for all s G [0, i] so that <£> y (s, y) = ^ ^(y) for s G [0, t] 
by (2.17). Since 

y) = 1 + tu x (y) + / (t - a) ^(s, y) . (2.23) 
J o 
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we find that in this case 



£„(*,!/) = l + tu x (y) + t -ul(y) = J (2 + tu x {y)) 2 . (2.24) 

In the remaining cases we have that u x (y) = — < — f for some t € (0,t\. Therefore (2.17) 
yields ip y (s,y) = \ u x {y) for s G [0,to) 5 while ip y (s,y) = for s 6 (ioj*]- From (2.23) we infer that 

e w (f,v) = 1-^ + ^=0. (2.25) 

The relations (2.24)-(2.25) confirm the monotonicity of the map y i-> £(t,y). Since £(0, x) = x, it 
remains to prove that linij^ioo y) = ±oo for any t > 0. Fix t > 0. Since u x € L 2 (IR), the 
Lebesgue measure of the set {y £ 1R : u x (y) < — j} is finite. On the complement C(t) of this 
set we obviously have u x (y) > — \ and thus £ y (t,y) > \ by taking into account (2.24). Therefore, 
given X2 > x 1} we infer that 



€(t,X2)-Z(t,x 1 )= [ Z y (t,y)dy> f £ v {t,y)dy>[ \ dy > 

Jx! J[xi,x 2 ]nC(t) J[xi,x 2 \ nc(t) 4 



X2 — Xl — Z(t) 



This proves the claim about the limiting behaviour of £(t, ■) at ±oo. While for times i up to the 
blow-up time T , given by (2.7), the map y i-> y) is an absolutely continuous diffeomorphism 
of the real line, for t > T this map is nondecreasing and onto but is not necessarily a bijection. 
Nevertheless, we would like to define the solution u by the formula (2.21) for all t > 0. 

To show that u is well-defined via (2.21), due to the monotone and surjective character of 
the map y i-> y), it is sufficient to show that if £(i,yi) = £(t,y 2 ) for some y 2 > yi, then the 
values of u given by (2.21) are also equal. Indeed, we must have that £(t,y) = £(i,yi) for all 
2/ S [2/1,2/2] and a glance at (2.24)-(2.25) confirms that u x (y) < — | for y 6 [2/1,2/2]- This means 
that for every fixed y 6 (2/1,2/2) we have ^(y) = — ^|^y for some io(y) G [0, i]. Consequently 
ip y (s,y) = \ u x (y) X[o,t (j/)]( s ) f°r s e [0)*] an d differentiation of the right-hand side of (2.21) yields 

/ /"* \ /•*o(y) 1 2 tnfty) 4 

% (fiw + 1 v) ds) = «.(,) + y o - ^ (»> ^-^g+f 

for ?/ € (2/1,2/2)- In particular, the values of the right-hand side of (2.21) are equal when evaluated 
at (i, 2/1) and at (t, 2/2) - This proves that tt is well-defined. 

The next step is to prove that for every t > 0, the map 2/ > 2/) is continuous on 1R with 
distributional derivative in L 2 (M). Given t > and 2/0 £ -K, let Iq = {x 6 iR : x) = 2/0}- The 
previously established properties of the map x 1— > £(i, x) ensure that Jo = [a, 6] for some a < b. For 
any sequence y n — > yo, choose x n £ M with x n ) = y n . If we show that min {|x n — o|, |x n — fe|} — ► 
as 71 — ► 00, by the continuous dependence on the y-variable of the right-hand side of (2.21), we 
infer that 

u(t, y n ) = u(t, £(t, n)) -> u(t, £(t, a)) = u(t, b)) = u(t, y ) 

since x n ) — ► £(t,a) = £(t,b) = yo- Thus y 1— > u(t,y) would be continuous at yo- If it would be 
possible that min{|x n)c — a\, \x nk — b\} > e > for a sequence — > 00, then 

\y nk - yo\ = \Z(t,x nh ) -£(t,a)\ = \£(t,x nh ) -£(t,b)\ > mm{y Q -(,(t,a-e), £(i,6 + e)} > 

must hold by the definition of [a, b] and the monotonicity property of the function x 1— > £(t, x). But 
this is a contradiction since y n — > yo as n — > 00. We therefore proved the continuity of the map 
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y i-> u(t, y) for every fixed y € M. Actually, a glance at the previous considerations confirms the 
continuity of the map u : M + x M —> M since £ : 1R+ x 1R — ► iR is continuous in view of (2.14). To 
show that for each t > the distributional derivative u x (t, •) belongs to L 2 (M), due to the absolute 
continuity of the nondecreasing surjective map £(£, •) : iR — > iR, we first show that at every point 
y = £(t,x) where £ x (t,x) > exists, u x (t,y) £ M exists. Indeed, at such a point y the right-hand 
side of (2.21), formally equal to u x (t,£(t,x)) ■ £ x (t,x), is differentiable with derivative 



;(x)+j <p y (t,x)ds = u x (x) + -u x (x), 



since in view of (2.24)-(2.25) we must have u x (x) > — | and ip y (s,x) = ^u x (x) for all s G [0,t]. 
Since £ y (t,x) = 1 + tu x (x) + ^- -u^(x), we infer that u x (t,y) exists, being given by the formula 

u x {t,y) = f = _ (2.26) 

l+te x (x) + ^4(x) l + I^a:) 

where y = £(£,«)• From (2.26) we deduce that for any interval [xi,^] where £ x (t,x) > a.e., we 
have 



u 2 x (t,y)dy= / u*(t,f(t,a;)) • ^(i,x)(ix 



r , Ux{x) , (i+tu x (x)+ t -uKx))dx= r 

J y> (i + \u x {x)) v 4 1 J y> 



u x (x)dx 



if Hi = £,(ti x i)i Vi = £(^£2) and if we take into account (2.24). Summing up over such intervals, 
we obtain that 

/ u 2 x (t,x)dx= / u 2 x (x)dx. (2.27) 

J M J {U x (x)>- 

In particular, the map t 1— > ||it x (i) •) Hi 2 (iR) i s nondecreasing on M + . Moreover, we also have that 



ul(t,x)dx= I 

-00 ^ {a 



r«*,v) 

■u x (x) dx 

{rre(-oo,{/]: it* (a) > - f } 



and 



00 



u x (t,x)dx = / u 2 (x)dx. 

£(t,y) J{xe[y,oo): u x (x)>- §} 



A comparison with (2.18) yields 

Furthermore, if ^(t, x) > exists, relation (2.20) ensures for y = £(t,x) the existence of 
u t (t,y), given by the formula 

Ut(t, y) = - u x (t, y) u(t, y) + v?(i, x) 
11 



obtained by differentiation and taking into acount (2.22). In combination with (2.28), this yields 

(u t + uux) (t, £(t, x ^ = \ J R si e n (^> x ) ~ c) C) d(, 

which is precisely (1.1) evaluated at (t,£(t,x)). In view of the previously established properties 
of the map x i-> £(t,x) we deduce that the constructed function u satisfies also the condition (iii) 
of Definition 2.1. Since the other properties required by Definition 2.1 were proved above, we 
conclude that u qualifies as a solution to (1.1) in the sense of Definition 2.1. This completes the 
proof of Theorem 1. <0> 



3 - A distance functional 

If u : M — > M is also bounded, in addition to being continuous and with distributional 
derivative u x 6 L (M), then the global solution u(t, •) constructed in Theorem 1 will be bounded 
at every fixed time t > 0. More precisely, in view of (2.19) and (2.21) we have that 



sup |ti(i,x)| < sup \u(x)\ + j / u 2 (x)dx. 

t>0,xeM x€M 4 J R 



Thus, if u : M — * M is a bounded continuous function with distributional derivative u x G L 2 (M), 
then at each fixed time t > 0, the solution u(t,-) to (1.1), constructed in Theorem 1, belongs to 
the Banach space X of bounded continuous functions / : M — ► M with distributional derivative 
f x G L 2 (M), endowed with the norm 



x =sup{\f(x)\}+([ fl{x)dx) 



The Banach space X seems suitable for (1.1) - see also [BZZ] where a construction similar to the 
one performed in Theorem 1 is presented. However, the map t *— > u(t, •) is generally not continuous 
from M + to X. Indeed, if for some r > we have that the set {x 6 M : u x (x) = — 7} is of 
positive Lebesgue measure, then a discontinuity occurs at time t = r for the map t i-> u(t, ■) € X 
since from (2.27) we infer that for t < r, 

/ u x (t,x)dx— / u^(r, x)dx> / tt^(x)cix>0. 

JiR JJR J{x€x: u x {x) = - ^} 

Our aim will be to construct a distance functional in the space of solutions to (1.1) with respect 
to which we will have both continuity with respect to time as well as continuity with respect to the 
initial data for the solutions to (1.1). More precisely, for non-smooth solutions the conservation 
law (2.2) is replaced by 

(v 2 ) t + (uv 2 ) x = -n, (3.1) 
where fx is the positive measure on the t-x plane defined as 



= If, ,1 u 2 x (y)dy 

'[{T{y), €(T(y),y))enj 
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for every open set 17 C M + x JR. Here T{y) is the blow-up time along the characteristic curve 
starting at y, namely 

T r y \ = / -2/u x (y) if u x (y) < , 

I oo otherwise. 

For any u £ X, we can use the semigroup notation S t u = u(t, ■) to denote the solution of (1.1) 
constructed in Section 2. Indeed 

S u = u, S t+S u = S t (S s u) . (3.2) 
To prove (3.2), we first show that 

Zi(t + s,y)=&(t,lii(s,y)), t,s>0,yeM, (3.3) 

where £2 is the characteristic built upon the initial data y i-> u(s,£i(s,y)). To check (3.3), we 
view both expressions as functions of t. At t = they are both equal to £i(s,y). For t > 0, 
differentiation of (3.3) yields 

ft+S rt 

u(y)+ ip 1 {r,y)dr = u(s,£ 1 {s,y))+ ip 2 (r,^(s,y)) dr (3.4) 
J ■/ 

in view of (2.20). We use (2.21) to express the right-hand side of (3.4) as 

«(!/)+/ <Pi(r,y)dr+ [ (p 2 (r,^(s,y)) dr. 
Jo Jo 

Therefore, to get (3.4), which yields (3.3) by integration, it suffices to show that 

/t+a rt 
ip 1 (r,y)dr = J ip 2 (r, £i(s, y)) dr. (3.5) 

To prove (3.5), we note that by (2.18), 

/ <P2(r,£ 1 (s,y))dr= - I [ sign(fi(s,y) - x) u 2 x (s,x) dxdr 

JO 4 Jo J{x: Ux (s,x)>-Z} v 7 

= 7/ / sign(^ 1 (s,y)-^ 1 (s,x))ul(s,^ 1 (s,x))dM(s,x)dxdr 

4 Jo J{x: u B («,€i(*.!r))>-2} V J \ J 

if we change variables x i-> £i(s,x). Taking now (2.9)-(2.10) into account, we infer that 

f <p 2 (r,^(s,y))dr = \ f [ sign(^ 1 (s,y)-^ 1 (s,x))u 2 x (x)dxdr 

JO 4 JO J{x: u B ( i ,€i(a,x))>-a} V 7 

sign(y — x) (x) dx dr 



1 '* 



4 ./0 J{x: u x (s,£, 1 {s,x))>-Z} 



since the function x i-> £i(s,x) is nondecreasing. But 

= 2 ^ X _ X > -- if and only if u x (x) > - 



2 + su x {x) r ' s + r 
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2y 

since the function y i-> is strictly increasing for y > — -, so that in the end we get 

2 + sy 

/ ¥>2(r,£i(s,2/))dr = - / / sign(y - x) u|(x) dxdr 

^0 4 JO ./{x: ^(x)>-^} 



i /• 

4 A 7 {2 



(3.6) 



sign(y — x) -u x (x) <ix dr 



A I I — o--\o ~/ ~"XV 

* is J{x: u x (x)>-^} 



where r = r + s. On the other hand, by (2.18), 

"t+s -1 rt+s 

! 

[x: u x (x)>-f} 



rt+s rt+s r 

\ tpi(r,y)dr = - \ / sign(y - x) u 2 x {x) dx dr 

Js 4 Js J{x: u x (x)>-2 : } 



so that (3.4) holds and (3.3) is proved. Knowing (3.3), to infer S t + S u = S t (S s u), it suffices to show 
that 

u(t + s, 6 (t + s, y)) = u(t, &(t, y))). 
But, by (2.21), the left-hand side is precisely 

rt+s rs rt+s rt+s 



dr, 



rt+s rs rt+s rt+s 

u(v)+ <p(r,y)dr = u(y) + ip 1 (r,y)dr + Pi(r,y) dr = u(s,^t(s,y)) + <fi(r,y) 

JO JO Js Js 

which, taking into account (3.5), equals to 

u(s,£t(s,y)) + f (P2(r,£ 1 (s,y))d8 = u(t,£2(t,£i(s,y))) 
Jo 

in view of (2.21). This completes the proof of (3.2). 

Notice that in general the map t t— > S t u is NOT continuous from [0, oof into X. It is thus 
interesting to identify some distance J(u, v) which is well adapted to the evolution generated by 
(1.1). More precisely, given an arbitrary constant M, in this section we shall construct a functional 
J(u, v) with the following property: For any initial data u, v £ X with 

\\u x \\i? < M , ||^x||l2 < M, 

the corresponding dissipative solutions u, v constructed in Theorem 1 satisfy 

J(u(t) , v(t)) < e CMt J(u, v). 

To begin the construction, consider the metric space 

X = (iR 2 x] -tt/2, tt/2]) U{oo} (3.7) 

with distance 

d(^(x, u, w), (x, u,w)^j = min | \x — x\ + \u — -£t| + Kq \ w — u)| , kq\tt /2 + w\ + kq\tt /2 + , 
d(^(x,u,w), oo^j = n \ir/2 + w\ . 

(3.8) 
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Here kq is a suitably large constant, whose precise value will be specified later. Notice that X is 
obtained from the metric space M 2 x [— tt/2, it/2] by identifying all points (x, u, — ir/2) into a 
single point, called "oo". 

Let M(X) be the space of all bounded Radon measures on X. To each function u G H} oc {lR) 
with u x G L 2 (M) we now associate the measure \i u G A4(X) defined as 

^({oo})=0, f(A)= I u 2 x (x)dx (3.9) 

^jxgiR: (x, tt(a;), arctan », (a;) ) 6 A | 

for every Borel set A C M 2 x } - tt/2 , tt/2] . 

As distance between two functions u, t> G X we now introduce a kind of Kantorovich distance 
J(u,v) related to an optimal transportation problem. Call T the family of all triples (-0,01,02), 
where 0i, 02 '■ JR l— * [0, 1] are simple Borel measurable maps (that is, their range is a finite number 
of points and the preimage of each such point is a Borel set) and ip : M t— > M is a nondecreasing 
absolute continuous surjective map. Assuming that 

4>!(x)ul(x) = ip\x) ■ 4> 2 {ip(x))vl(i/j(x)) fora.e. xeM, (3.10) 

we define 

arctan u x (x)) , (ip(x), v(ip(x)), arctan v x (ip(x))*J ■ 0i (x) u 2 (x) dx 
+ /rf((x, «(*),«■*»«.(,)), »).(1 -#,(.)) 

+ .(««)), , oo) • (1 - * 2WW ))„ ;W x))^(x )(J x. 

(3.11) 

Observe that (-0,01,02) can be regarded as a transportation plan, in order to transport the 
measure fj, u onto the measure \i v . Since these two positive measures need not have the same total 
mass, we allow some of the mass to be transferred to the point oo. More precisely, the mass 
transferred is (1 — 0i) ■ [i u and (1 — 02) • /i l '. The last two integrals in (3.11) account for the 
additional cost of this transportation. Integrating (3.10) over the real line, one finds 



<fa.(x) u x (x) dx = / (p 2 (y)v x (y)dy . 

M JR 

We can thus transport the measure 0i^ u onto (f>2 /J- v by a map ^ : (x, u(x) arctan -u^x)) i-> 
(y, v(y), arctan v x (y)), with y = ip(x). The associated cost is given by the first integral in (3.11). 
In this case the measure 02 ^ v is obtained as the push-forward of the measure 0i// u . We recall 
that the push-forward of a measure \i by a mapping \t is defined as (^#/jl)(A) = fi(^~ 1 (A)) for 
every measurable set A. Here ^~ 1 (A) = {z : SSf(z) G A]. 

We now define our distance functional by optimizing over all transportation plans, namely 

J(u,v)= inf {jW'+^fav)} (3.12) 

where the infimum is taken over all triples (-0, 0i, 02 ) G ^ such that (3.10) holds. 

To check that (3.12) actually defines a distance, let u, v , u> G X be given functions. 
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1. Let us show that J(u,v) = J(v,u). In order to do this, it is enough to prove that for every 
triple (ip, (pi, <p2) G J 7 satisfying (3.10) and every e > 0, there is a triple (rj, tpi, ip%) G T satisfying 
(3.10) such that rj : M — > IR is a strictly increasing absolutely continuous bijection and 

jfo.^ 2 )( U)U )_ i 7(iMi.fc)( Ujt; ) < e . (3.13) 
Indeed, given (ip,<pi,(p2) G T satisfying (3.10), define ip = (p\ = <f2, <p2 = '■Pi- The properties 



u,v) 



of 77 ensure the absolute continuity of tp (see [N]) so that we obtain ji^^ 1 ^ 2 ) (v, u) = J^^ 1 ^ 2 ) 
by performing the change of variables x 1— > 77(2;). Since e > was arbitrary, we infer that J(v, u) < 
J(u,v). Interchanging the roles of u and v we get J(u,v) = J(v,u). 

To prove (3.13), it is convenient to view ip : IR — > M as a maximal monotone multifunction 
ip : M i-> V(M) with domain and range IR. Here V(M) is the family of all subsets of JR. The 
conditions for a multifunction F : IR 1— > V(M) to be maximal monotone with domain and range 1R 
may be explicited as follows [Z]: 

- for every x G IR, the set F{x) C M is nonempty (i.e. the domain of F is M); 

- for every y £ IR there is at least some x G IR with y € F(x), expressing the fact that the 
range of F is M; 

- there are no couples (xi,yi) and (x2,y2) with y\ G F(xi) and y 2 G F(x 2 ) such that xi < x 2 
and 2/2 < Di, meaning that F is monotone); 

- if we associate to F its graph {(x, y) G M 2 : y G F(x)}, then this graph has no proper exten- 
sion satisfying the first three properties (condition defining the maximal monotonicity property) . 

We recall some important features presented by such maps [A A, Z]: 

- the set F(x) is an interval of the form [a x , b x ] with a x < b x for all x G JR and a x = b x for 
all x G M, except perhaps an at most countable set (so F is singlevalued with the exception of at 
most countably many points); 

- F is a.e. differentiable, that is, for almost all x G JR there exists F'(x ) G IR such that 

lim y - F(x ) - (x - x ) F'(x ) = Q _ 

x^x a ,yeF{x) X — X 

- we can define the inverse F _1 : IR — > V(M) of F by asking y G F _1 (x) if and only if x G F(x) 
and F _1 is again a maximal monotone multifunction with domain and range M. 

Since the multifunction ip -1 is maximal montone, let {y n } be the (at most countable) set of points 
where it is multivalued, that is, ip~ x {y) = [a n ,b n ] with b n > a n . Then ip(x) = y n for x G [a n) ^n] 
and, ip being absolutely continuous, tp x > a.e. on M — \J n [a n ,bn] since ip is strictly increasing 
on this set. Given 7 > 0, the absolute continuity of ip : M — > M allows us to choose some S > 
such that the total variation of ip over the union of disjoint closed intervals with the sum of their 
lengths less than 5 is less than 7 cf. [BGH]. On each interval [a n — <pr, b n + ^-] we replace ip with 
the linear function r\ which takes the values ip(a n — ^-), respectively ip(b n + ^-) at the endpoints. 
By the way a n and b n were defined, we know that ip(b n + ^r) > ip{a n — ^-) so that r/'(x) is a 
positive constant on [a n — Jr, b n + ^] with 



/ r/(x) dx < [Wn + —) - ip(a n - — )J < 



E 

n " " n 2" 



7- 



Setting r/(x) = ip(x) for x [a n — 6 n + Jr], we obtain a strictly increasing absolutely continuous 
bijection 7] : H ^ M. Let us now show that the triple (ry, ¥>i,<£>2) G J 7 satisfies both (3.10) and 
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(3.13), where tpi, f 2 are defined by setting <pi(x) = for x G [a n — b n + ^] and (fi(x) = (f>i(x) 
for x g [a n - Jr, b n + while <p 2 (r](x)) = 4> 2 (ip(x))) for x g [a n - b n + and <£ 2 (f?(z)) = 
for x G [a n — b n + ^-]. On the complement of the set Unt a ™ ~ ^jb n + ^] relation (3.10) clearly 
holds a.e. for (r/,(pi,(p 2 ), being unmodified from (3.10) for (ip,(j)i,(f) 2 ). If x € [a n — ^r,6 n + ^r], 
then (3.10) for (r],ipi,ip 2 ) holds again since both sides are zero as ipi(x) = (p 2 (r)(x)) = in this 
case. Finally, to check (3.13), notice that if we denote 

E 5 = U{K - — , a n ] U [b n ,b n + — ]}, A = [j[a n ,b n ], 

n n 

then 

J( w *)(m^)-J< Wi *'(ii,d) <2k i f u 2 x dx + 2h tt [ v 2 x (r](x)) rf(x) dx. (3.14) 

Je 5 Je 5 u a 

Indeed, the distance d is less that 2k 7t and the integrands in J^'^ 1 ^ 2 ^ (u, f ) and 
agree on the complement of the set UrJ a ™ ~~ ^T>b n + ^-] by definition. Also, for a.e. x € [a n ,6 n ] 
we have 4>i(x)u 2 = by (3.10) as ip'(x) = 0, and <pi(x) = by its definition. We obtain (3.14). 
Since the absolutely continuous map r\ maps E$ U A into a set of Lebesgue measure less than 7, 
and u 2 , v 2 € L x (iR), from (3.14) we infer (3.13) by choosing 5 > and 7 > small enough. This 
completes the argumentation needed to show that J(u,v) = J(v,u). 

2. Choosing ip(x) = x, (f>i(x) = <j) 2 (x) = 1, we immediately see that J(u,u) = 0. Moreover, we 
have J(u, v) > if u 7^ v. To check this, note that J(u, v) = implies that there is a sequence 
(V> n ,0i, 4%) along which J^"^" v) -» 0. The second term in (3.11) yields 

(- + arctanu a; (x) N ) (l-^(x))u 2 x (x) ^0 in L l (M), 



.2 

so that along a subsequence (1 — (j)i k )u 2 — ► a.e. on iR since > —00 a.e. On the set 
5 = {x 6 iR : u x (x) 7^ 0} we therefore have </>™ fc — ► 1 a.e. Moreover, the first term in (3.11) forces 

</>™ fc (x)u 2 x (x) ■ min j|x - ip nk (x)\ + \u(x) - v(ip nk (x))\ + k | arctanu x (x) - arctanu a; (^ nfc (x))|, 

k ^ + arctanu^x) + ^ + arctanu x (V' n ' c (x))J | — > in L x (iR). (3.15) 
Since u x > — 00 a.e. ensures 

7T 7T 7T 

— + arctanti x (x) + — + arctan (V> n,c (x)) > — + arctanu x (x) > a.e. on M, 

we infer from (3.15), by passing to another subsequence, that 

|x - ip Uk (x)\ + |u(x) - v(tp nk (x))\ -> a.e. on S 1 . 

In view of the continuity of v, ^™ fc (x) — > x a.e. on 5 guarantees -u(^ nfc (x)) — > -u(x) a.e. on S so that 
u = v a.e. on S since also v(ip nk (x)) — > u(x) a.e. on S 1 . Repeating this argument with the roles 
of u and w reversed, we find that u = v a.e. on the set {x £ M : v x ^ 0}. Combining this with 
the previous conclusion, we have u = v a.e. on the complement of the set {x € M : u x = v x = 0}. 
Since u x , v x € L 2 (M), this is possible only if u = v on M. Thus J(it, v) = if and only if u = v. 
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3. Finally, to prove the triangle inequality, it suffices to show that for every choice of (0 , 0i,0 2 ) 
satisfying (3.10), and of (ip$ , 4 , 2 ) satisfying (3.10) for (v,w), the triplet (tp, 0i, 2 ) defined by 

00*0 = V s (^ b (x)), <hix) = 4>\{x) ■ <l>W{x)), My) = 4(y) ■ <t>W(?)), 

satisfies (3.10) for (u,w) and 

jW-*»^)(u,«)) < J^'^\u,v) + J^^'^)( W)U; ). (3.16) 

Notice that composing the relation (3.10) for (v,w) a.e. to the right with ip^ , and multiplying 
the outcome by cj>\ o 0^ ■ (-0 b )', we infer that (3.10) holds a.e. on M for (u,w) with our choice of 
(0, 0i, 2 ) and we can now concentrate on proving (3.16). 

To simplify matters, we introduce the following notation 

Pi = (x, u, arctantt x ), P 2 = (0 b , t> o b , arctanfa, o b ), P 3 = (ip, w o ip, arctanu> x o ip), 

mi=u 2 x , m 2 = v 2 x o 0; b ■ (0 b )', m 3 = wl o ij) ■ ip' . 
The relations of type (3.10) yield then that a.e. on M, 

(f>\ • mi = 4> 2 ° ^ ' m 2, 4>\ ° "4 • m 2 = $2° ip ■ m 3 , 0i • mi = 2 ° V> • m%. (3.17) 

Also, 

jW'+^lu, w) = J {d(Pi, P 3 ) • ^imi + d(Pi, oo) • (1 - 0i)mi + d(P 3 , oo) • (1 - 2 o 0;)m 3 } dx, 

jW» k .*i.^)( UjU ) = y |d(Pi,P 2 ) • b mi + d(Pi, oo) ■ (1 - b )mi + d(P 2 , oo) • (1 - 2 o 0; b )m 2 } dor, 

jW»'.*"i = /" {d(P 2 ,P3)-0 S 1 o^ b m 1 +d(P 2 ,oo)-(l-0 S 1 o^)m 2 +d(P3,oo)-(l-0»o^)m 3 }dx, 

the last relation being obtained after the change of variables x i-> 0; b (x) in the integral. We will 
prove (3.16) by deriving an appropriate inequality valid a.e. pointwise between the integrands in 
the previous expressions. Since 

(l-0 2 o^)(l-0»o^)>O, 

we have 

1 - (j) 2 o 0; b + 1 - 0} o V b > 2 o V b (l - 0? o b ) + A o 0^ b (l - 02 o 0; b ). 
Multiplication of both sides by d(P 2 , oo) • m 2 leads to 

d(P 2 , oo) • (1 - b o 0, b )m 2 + d(P 2 , oo) • (1 - 4 o V b )m 2 > d(P 2 , oo) • b (l - 0» o ^)mi 

+ d(P 2 , oo) • 4 ° ^(1 " 2 ° (3-18) 
in view of (3.17). Multiply now the inequalities 

d(Pi,P 2 ) -d(Pi,oo) + d(P 2 ,oo) > 0, d(P 2 ,P 3 ) -d(P 3 ,oo) + d(P 2 ,oo) > 0, 
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by 4>\(1 — 4>\o ip b )m 1 , respectively <f>\ o 0; b (l — 4 2 o ■0 b )m 2 , and add them up. The outcome yields 
in combination with (3.18) that 

d(p 1 ,p 2 ) • - 4o^) mi - d(Pi,oo) • ^(i - 0? o vVi + ^(p 2 ,p 3 ) • 0! o v b (i - 2 o vV 2 

+d(P 2 , oo) ■ (1 - 4 o 0, V2 + d(P 2 , oo) ■ (1 - 4 o ^)m 2 > d(P 3 , oo) • <j>\ o/(l-^o ^)m 2 . 
Adding to both sides the quantity 

d(P u oo) • mi + d(P 3 , oo) • m 3 + d(P u P 2 ) ■ <j>\ ■ <f>[ o^-mi- d(P x , oo) • </>{ ■ 0? ° ^ ■ ™i 

+d(P 2 , P 3 ) • 4i -4 2 o^ .m 2 - d(P 3 , oo) • 4 2 o i\> ■ m 3 

we deduce by (3.17) that the integrand of J^'^'^(u,v) + w), equal a.e. precisely 

to the left-hand side of the new inequality, is a.e. pointwise larger than 

d(P 3 , oo) • 4i o V b (l - 02 ° ^ V2 + d(P u 00) • mi + d(P 3 , 00) • m 3 + d(P u P 2 ) ■ 4 ■ 4i ■m 1 

-d(Pi,oo) • b ! • 0« o 0, b • mi + d(P 2 ,P 3 ) • 4i o^-^o^. m2 - d(P 3 ,oc) • 2 o V • m 3 . 

Taking into account (3.17) and the definition 0i = 4 ' 4\ we see that the above expression 
equals 

d(P 1 , 00) • (1 - 0i)mi + d(P 3 , 00) • (1 - 02 o V)m 3 + (d(P u P 2 ) + ^(P 2 , P 3 ) • 0imi) 

> d(Pi, 00) • (1 - 0i)mi + d(P 3 , 00) • (1 - 02 o V)m 3 + d(P 1: P 3 ) ■ farm. 

The lower estimate is a.e. precisely the integrand in jW'* 1 '* 2 ' (u, u;) and (3.16) holds. The proof 
that J satisfies the triangle inequality is therefore completed. 

In the remainder of this section we examine how the distance </(•,•) behaves in connection 
with solutions of the equation (1.1). 

Continuity w.r.t. time. Let t i-> u(t) be the solution of (1.1) constructed in Section 2. For any 
fixed t > 0, we define a transportation plan of fi u to fi u ^ by setting 

#r) 0i(x) = |j !| T(x)<t! ^(x) = l. (3.19) 

Relation (3.6) follows from (2.9)-(2.10) on {T(x) > t} and from (2.25) on {T(x) < t}. The cost of 
this plan is estimated by 

J (w "* 2) (m, u(t)) < [ \\x-£(t,x)\ + \u(x)-u(t,Z(t,x))\ 

J{T{x)>t] L 

+ K \ arctan -u x (x) — arctanUa;(t,^(t,x))|| u x (x) dx (3.20) 
+ / |"7r/2 + arctan-Ua;(x)| u x (x) dx . 

J{T{x)<t} 
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By (2.4) we have that a.e. 

d_ 

dl 



arctan (t, £(t,x)) 



f t u x (t^(t,x)) 



l + u*{t,Z(t,x)) 



< 



(3.21) 



An integration on [0, t] yields 

| arctan-u x (x) — arctanu x (t, £(i,x))| < ^, t > 0. (3.22) 
On the other hand, using (2.20), we get 

\x-£(t,x)\ < t\u(x)\+ f (t-s)\<p(s,x)\ds <t\u(x)\ + t — [ u 2 x (x)dx, t > 0, x G M, (3.23) 

Jo ° Jm 

if we take into account (2.19). From (2.21) and (2.19), we also infer 

\u(x)-u(t,Z(t,x))\ < f \<p(s,y)\dy<l I u 2 x (x)dx, t>0,x€M. (3.24) 

Jo 4 Jr 

To estimate the last term in (3.20), notice that 

2 2 
{x G M : T(x) < t} = {x G M : - — - - < t} = {x G M : uJx) < - -}, t > 0. (3.25) 

u x (x) t 

TT 

Furthermore, since lim x{ — h arctanx) = —1, there is a constant c > such that 

X — > — oo 2 

TT C 

< - + arctany < — , y < -1, 

2 |y| 



whereas 



so that 



7T 



+ arctan y 



y 2 < tt if - 1 < y < 0, 



7T 



— + arctan -u^x) ^(x) < 7r + c |-Ua;(x)| if ^(x) < . 



On the other hand, if u x {x) < — - , then t 2 u 2 (x) > 4 so that 



{T(x)<t} 



ldx < — 
~ 4 



{T(x)<t} 



u x (x) dx. 



From (3.25)-(3.27) we infer that 



f n I \ -2 / \ 7r * 2 n - n2 /" 

/ h arctan -u^x) u x (x) dx < —— \\u x \\ L2 + c I 

J{T(x)<t\ 2 4 Jr 



'{T(a;)<t} 

< \\^x\\h 4 - ( 



{T(x)<t} 



\u x (x) \ dx 



([ ldx) 2 (f u 2 x {x)dxY < 7 ^— 

J{T(x)<t} J ^J{T(x)<t\ 1 



(3.26) 



(3.27) 



>{T(x)<t} 
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U x\\L2 + — \\U X \\ 2 L 2- 



By (3.20), (3.22)-(3.25) and the previous inequality we conclude 

J(^^)(fi, u(t)) < + + ||u|| L - + ||^||| 2 ) 1 1|^||| 2; t > 0. (3.28) 

It is now clear that each semigroup trajectory t t— > 5 t -u is Lipschitz continuous as a map from 
[0,oo [ into the metric space X equipped with our distance functional J. The Lipschitz constant 
remains uniformly bounded as u ranges over bounded subsets of X. 

Continuity w.r.t. the initial data. We now consider two distinct solutions and study how the 
distance J(u(t) , u(t)) varies in time. Recall that the solution u = u(t,x) is computed by (2.20)- 
(2.22), also in the case where the gradient blows up. The same formula of course holds for u. Let 
(^0)01,0)02,0) be an optimal transportation plan of the measure /U u (°) to the measure fi u( -°\ In 
view of the approximation property established in (3.13), we can restrict our attention to the case 
when -00 is strictly increasing on M. For any t > 0, we define a transportation plan (-0*, 4>\, 2 ) of 
the measure fi u ^ to fi u ^ as follows: 

^{^t,y))=i(t, y) for y = My), 



0i(^y)) = 




if T(y)>t and T(y) > t for y = o (y), 
otherwise, 



0* y)) = (02,o(y) if T(y)>t and T(y) > t for y = O 1 (y) , 

[O otherwise. 

If initially the point y is mapped toy = ipo(y), then at any later time t > the point y) along the 
it-characteristic starting from y is sent to the point £(t, y) along the ^-characteristic starting from 

y = ipo(y)- We thus transport the mass from the point (^(t,y) , u(t, £(i,y)) , arctann^ (t, £(t, y))^ 

to the corresponding point (ji(t,y) , u(t, £(i, y)) , arctan{t x (f,£(i,y))^ with y = ipo(y), ex- 
cept in the case where blow up has occurred within time t along one (or both) characteristics 
£(-,y), C(-,y). In this later case, the mass is transported to the point oo. 
To check (3.10), it suffices to show that a.e. 

Since the relations y = i/j (y), ^*(£(t,y)) = £(t, y), and (V>*)'(£(t, y)) • ^(t, y) = | x (t, ipo(y)) ■ il>' (y) 
all hold a.e., the desired identity holds a.e. on the complement of the set {y : y = i/>o(y), T(y) > 
i, T(y) > i} where both sides equal zero since 0i(£(i, y)) = 02(£(*>y)) = 0. The identity holds 
also a.e. on the set {y : y = ipo(y), T(y) > i, T(y) > t} since there, in view of (2.9)-(2.10), it 
practically amounts to relation (3.10) for (01,0) 02,0) Vto)- 

In the following, our main goal is to provide an estimate on the time derivative of the function 

J(t) = jW*>+*>+*\u{t),u{t)). 

Throughout the remainder of this section, by {T(y)<i} we understand the set of all y 6 M such 
that MV) = {y} and f(y)>t. Since u 2 x (t,^t,y)) ■ £~(t,y) = u 2 x (y) on {T(y) > t} by (2.9)-(2.10) 
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and u 2 x (t, i(t,y)) ■ = u 2 x (0,y) on {T(y) > t}, while Y>*(£(i,y)) = i(t,y) for Mv) = iv}, 

performing the change of variables y i— > £(i, y), we see that 

J(t)= [ min \m,y)-i(t,y)\ + \u(t^(t,y))-u(t,i(t,y))\ 

J{T(y)>t,T(y)>t} 

+ k | axct&nu x (t,£(t,y)) - axctaa.u x (t,£(t,y))\, 

k (n + arctan u x (t, £(t, y)) + arctan u x (t,£(t,y)f)} (j> lj0 (y) u 2 x (y) dy 

+ k / + arctan u x (t,£(t,y))) [1 - 4>i,o(y)) u 2 x (y) dy 

J{T(y)>t,f(y)>t} V 2 ' V ' 

+ k / (77 + arctan« a .(t,^(t,y))J u 2 x (y)dy 

J{T{y)<t or T(y)<t} v ^ 7 

+ k / ( J + arctan^ y))) (l - ^ 2 ,o(y)) u 2 x (0,y)ijj' (y)dy 

J{T(y)>t,f(y)>t} V 2 ' V ' 

+ k + arctan tt x (t, y))) u 2 x {0,y) ip' (y)dy. 

HT(y)<t or T(y)<t} v ^ 7 

To simplify notation, let 

S(t) = {T(y)>t,f(y)>t}, S c (t)=M-S(t), (3.29) 

= min \\u(t,£(t,y)) - u(t,i(t,y))\ + k \ arctanu x (t,£(t,y)) - arctan u x (t, y))| 

- k o(tt + arctanu a; (t,^(t,y)) + arctan {^(i, £(i, y))) |. (3.30) 
Since for h > we have 

S(t + ft)cS(f), S c (i)cS c (f + li), (3.31) 

we deduce that 

J(t + h)- J(t) = [ (E(t + h,y)- E(t, y)) faM u 2 x (y) dy 
Js(t) v 7 

-/ E(t + h,y)(f> 1 , (y)ul(y)dy 

JS(t)\S(t+h) 
/ + *,«* + *.»)) —tan 

/ + + —tan (t,|(t,jl))) (i-0 2 ,o(y)) ■4(0, y) ^o(y)dy 

-«o / (77 + axctan« a .(t + /i,^(i + /i,y))j (l - 0i,o(y)) u 2 x (y)dy 

JS(t)\S(t+h) KZ / \ / 

-k / ( J + arctan u x (t + /i,c|(t + /i,y))) ( 1 - </> 2 ,o(y) ) ^(0, y) ^o(y) d 2/ 

JS(t)\S(t+h) v/ / \ / 

+K0 / (arctan„,( t + k ,« t + k,y)) - arctan ,))) 



'S=(t) 
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+k / (arctanu x (t + h,£(t + h,y))-arctanu x (t,£(t,y))) ul(0,y)ij)' o (y)dy 
Js c (t) v 7 

+k / (77 + aTct&nu x (t + h,£(t + h,y))) u 2 x (y)dy 

J S^{t+h)\S^{t) K/ 7 

+k I (^ + aTctanu x (t + h,i(t + h,y)))ul(0,y)^ o (y)dy. (3.32) 

JSc(t+h)\Sc(t) yi 7 

Noticing that S(t) \ S(t + h) = S(t) n S c (t + h) = S c (t + h)\ S c (t), we see that the combination 
of the fifth and ninth terms above, with that of the sixth and tenth, added to the second term, 
amount to 

k j + arctan u x (t + h,£(t + h,y)f) ^(y) u 2 x (y) dy 



S(t)\S(t+h) 



+k / + arctanu x (t + h, £(t + h,y))) (f> 2 ,o(y) ul(0,y)^' o (y) dy 

JS{t)\S(t+h) yi 7 

- I E(t + h,y)(f) 1: o{y)u 2 x (y)dy 

JS(t)\S(t+h) 

= k q (tt + arctanu^i + h,£(t + h,y)) + arctanu^t + h,£(t + h,y))) (/>i,o(y) u 2 x {y) dy 

JS{t)\S{t+h) v 7 

-[ E(t + h,y)^ (y)u 2 x (y)dy (3.33) 

JS(t)\S(t+h) 

by (3.10) for (0i, o , 02,o, ifo). In view of (3.29)-(3.30), on S(t) \ S(t + h) we have that 

E{t + h,y) = 7T + arctanua,(t + h,£(t + h, y)) + arctanu a; (t + h,£(t + h,y)) 

since at least one of the expressions u x (t + h,£(t + h, y)) and u x (t + h,£(t + h, y)) is precisely — oo 
on this set. Thus the whole expression (3.33) is identically zero. Therefore (3.32) yields 

j{t + h )- J(t) = [ (E(t + h,y)- E(t, yj) fa fi (y) u 2 x (y) dy 
Js(t) v 7 

u x (t,t(t,y))) (l-^i,o(y)) u 2 x (y)dy 
+ Ko J s (arctanu x (t + h,£(t + h,y))-arctanu x (t,£(t,y))) (l - <h,o(v)) 

Js c (t) v 7 

+k / ( arctanu^ + + h,y)) - axctanu x (t,£(t,y))) ul(0,y)ip' o (y) dy. (3.34) 
Js c (t) v 7 

In view of (3.29), we have S c {t) = {T(y) < t} U {f(y) < t}. On the set {T(y) < t} we have 
arctanw x (i + h,£(t + h, y)) = arctan u^i, £(i, y)) = — oo so that in the fourth term in (3.34) only 
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the integral over {T(y) > t, T{y) < t) might have a nonzero contribution. Thus the second and 
fourth terms in (3.34) combine to 

k q / ( a,rctanu x (t + h,£(t + h,y))-arcta,nu x (t,£(t,y))) ( 1 - 4>i, (y) ) u 2 x (y) dy 

J{T(y)>t} v J \ J 

+k / (avctanu x (t + h,C(t + h,y))-arctanu x (t,^(t,y)))(f) 1: o(y)ul(y)dy. (3.35) 

J {f(y)<t<T(y)} V ' 

Similarly, the third and fifth terms combine to 

i / ( aictanu x (t + h,£(t + h,y)) - aictanu x (t,£(t, y))) ( 1 - (f> 2 ,o(y) ) ^(0, y) i/)' (y) dy 

J{f(y)>t} V / \ / 



Kq 



+k I ( arctan u x (t + h,£(t + h,y)) - arctan u^i, £(t, y)) J <fo,o(y) u x (0,y) ip (y) dy. 

J{T(y)<t<f{y)} V J 

(3.36) 

To transform suitably the first term in (3.34), let us denote by E x {t,y) the first expression in the 
minimum (3.30), and by E 2 (t,y) the second. If E(t,y) = E 2 (t,y), then 

E(t + h,y) — E(t, y) < n ^ arctan -u^i + h,£(t + h, y)) — arctan u x (t, £(t, y)) 



+ arctan u x (t + h,£(t + h,y)) — arctan u x (t, £(t, y))^j , 



(3.37) 



since E(t+h, y) < E 2 (t+h,y). On the other hand, if E(t, y) = E 1 ^, y), then the triangle inequality 
and the relation E(t + h,y) < E 1 (t + h, y) ensure that 



E(t + h,y)- E(t, y) < £(t + h,y)- £(t, y) + £(t + h,y)- |(t, y) 



+ 



u(t + h, £(t + h,y))- u(t + h, i(t + h,y))- u(t, £(t, y)) + u(t, £(t, y)) 



(3.38) 



arctan u x (t+h, £(t+h, y))— arctan u x (t, y))+arctanu x (t+h, £(t+h, y))— arctan u x (t, £(t, y)) 



Letting h J, in (3.34), and taking into account (3.38) and the considerations preceding it, we 
deduce that 



lim sup 

h[0 



J(t + h)- J(t) 
h 



f d d 

< Ko J (t)+ / c/>i,o(y)ul(y) jA&y)- jJ&My)) 

J S(t~) 

dy 



dy 



d 



d 



+ J My) u 2 x (y) j t u(t, at, y)) - j t mt- > •<> (//)) ) 



+'■-■() / <f>i,o(y) u l(v) ^ arctan 2/)) - ^ arctan (4 (/. (',)(/;))) 



S(t) 



dt 



dy (3.39) 



where 



'[T(y)>t] V J 



— arctan u x (t,£(i, 2/)) u x (y) dy 
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+ 



J[T{y)>t] v 



2. o(//)) j t arctan u x (t,£(t,y)) u 2 x (0,y)dy 



+ 



I, 



[T(y)<t<T(y)] 



fafliv) [4 arctan u x (t,£(t,y)) \ u 2 x (y) dy 



-dt 



+ 



[d 



T{y)<t<T{y)] 



4>2,o(y) arctan u x (t, £(t, y))J u x (0, y) dy < 0, 



the last inequality being true by (2.4). 

Before proceeding with the further analysis of (3.39), we establish a few a priori bounds. From 
(2.22) we get 

4 y)~i f (*> y) = «(*> y)) - £)) • (3.40) 



Also, note that if w = arctan z(t) and z = — — , then 



2 + 2z 2 



- sin v. 
2 



Since | sin 2 a — sin 2 /3| < \a — (3\ by the mean-value theorem as |(sin 2 z)'| = 2|sinz cosz| 
| sin (2z)\ < 1, we infer three useful facts if we set z = u x (t,£(t, x)). First of all, 



arctanu x (t,£(t,y)), arctan u x (t, u(t, £(t, y)) < 0. 



(3-41) 



Secondly, 



arctan u(t, £(t, y)) — arctan u(i, £(t, y)) < ^ arctan (i, £(t, y)) — arctan« x (t, y)) 



7T 1 

Furthermore, if arctanz < , then sin(arctanz) G [—1, — 1, so that 

~ 4 v^ J 



(3.42) 



arctan u x (t,£(t,y)) < - j if arctan -u x (i, £(t, y)) < - ^. 



(3.43) 



On the other hand, using first (2.21) and then (2.18), we have 



^u(t,£(t,y )) - ^u(t,i(t,y )) = <p(t,yo)-<p(t,y )) 
/ sign(y - y) u 2 x {y) dy - / sign(y - y) u 2 x (0, y) 

JiT(v)>t\ J{f(y)>t} 

yo 



'{T(y)>t} 

i /•»(■ 
4 



/s/o ryo 
«*(y)X[T(y)>t]dl/- / u 2 x (0,y)x [ f {y)> t]dy 
-OO J —DC 

/>00 /"OO 

/ «x(i/)X[T(y)>t]dy- / «*(0,lf)X[f(fi 
Vn J Vn 



T(y)>t] 



Uo 



u x (y)x[T( y )>t]dy - / fix(o,Vo(y))X[f(s)>t] ^o(y) d v 
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P CO f'OO 

I u 2 x (y)x[T{ y )>t]dy - / ^(0,V'o(?/))X[T(y)>t]^o(y) d y 
Jya -'vq 

1 f ya / \ 

' ,o(v) + ^,o(y))ul(y)x[T( y )>t]dy 

i - foAMv)) + 4>2,o(My)))u 2 x (o,My)) X[f(y)>t] ^0(2/) d y 



4 

yo 



Va 

1 



- / (i - 4>i,o{y) + 0i,o(y))^(y)x[T( J/ )>t]^y 
- / (i-(f>2AMy)) + hAMy)))^l{o,My))x [ f(y) > t]^o(y) d y 

Jyo 

after performing in the next to the last step in two of the integrals the change of variables y = ipo(y). 
Since (3.10) for (^0,^1,0,^2,0) ensures that <£i, (y) u 2 x (y) = (fo,oCPo(y)) m£(0, V>o(y)) ip' (y) a.e. on 
5(f), we deduce that 

— u(t, f (t, yo)) - — u(t, i(t, y )) 



1 

< - 

~ 4 



dt v '" uy/ di 
1 - My)) u 2 x (y) dy+- - [ (l - 02, o (y))^(O, y) dy 

[T(V)>*] v 7 4 ^[f (v)>t] V 7 

+7/ <t>iAy)v 2 x (y)dy+] [ cj) 2 Ay)ul{^y)dy. (3.44) 

4 ■/[T(v)<t<T(»)] 4 ^[T(j,)<t<f(y)] 

Let us now introduce the following sets 

7T 7T 

Si = {y € 1R : arctan u x (£, y)) < — — and arctan ^(i, £(i, y)) < — — }, 

7T 7T 

£2 = {y € iR : arctan u^i, ^(i, y)) > — — and arctan u x (t, y)) > — 

7T 

5 3 = {y£M: arctan u x (t,£(t,y)) > --> arctan tt x (t,£(i, £))}, 

7T 

54 = {y € M : arctan u x {t,^(t,y)) > — — > arctan u x (t, £(i, y))}- 

The integral on the right-hand side of (3.44) over Si is, in view of (3.43), bounded from above by 
l^o(*)| = — Jo(t). The integral over S2 is, in view of the formula for J(t) preceding relation (2.29), 

bounded from above by ^Q-. To evaluate the contribution over the integral over S3, notice that 

7TK 

the same formula for J(t) yields 



J{t)>K / ( arctan u x (t,£(t,y)) - arctan u x (t,£(t,y)) ) fa^y) u 2 x {y) dy 

Js(t)nS 3 v 7 

+ K o / (? + arctan u x (t,£(t,y))) (l - 0i,o(y)) u 2 x (y)dy 

Js(t)ns 3 v/ / \ / 

+ k (77 + &rcta.nu x (t,£(t,y))) u 2 x (y) dy 

is c (i)ns 3 v 1 7 

+ K ° / ( J+arctanu^(t,|(t, y))) (l- ( /) 2i o(y))^(0,y)V'o(?/)^y 

Js(t)ns 3 v/ 7 v 7 

+ k / (77 + arctan ^(t, £(t, y))J u|(0,y) i/)' (y)dy. 

JS c (t)nS 3 Kl 7 



26 



Using (3.10), the sum of the first term and the fourth term is larger than 

JS(t)ClS 3 v/ 7 

f ul(0,y)^(y)dy=^f f u 2 x (0,y)dy 

JS(t)nS 3 4 Js(t)nS 3 



> Kq7T 



>T / (i-My))ul(o,y)dy = ?f- [ (i-<hAvj)*lM)dv, 

4 JS{t)ns 3 v 7 4 i[f(i)>i]ns 3 v 7 

with the last equality enforced by 63 C [T(y) > t]. The second term is bounded from below by 



kqtt 



/ (1 - 4>i,o(y)) u l(y)dy = / (i-0i,o(2/))*4(y)<% 

JS(t)nS 3 v 7 4 J[T(j)>t]nS 3 v 7 



4 Js(t)nS 3 V 7 4 i[T( s )x]nS 3 

■^T / f 1 - faAvi) u 2 x (y) dy= ^ / (l-0i,o(y))^(y)^y 

k vt 



/ u 2 x (y)dy+-^- \ (f> h0 (y) ul(y) dy 

J\T(y)>t>f(y)]nS* 4 J\T(y)>t>f(y)]nS 3 



'[T(y)>t>T(y)]nS 3 * J[T(y)>t>T(y)]r\S 3 

if we recall (3.10). The third term is bounded from below by 



/ *Uv)dv>*£[ ul{y)dy. 



4 Js-(t)ns 3 " 4 V[T(j/)>t>f(j/)]ns 3 

Summing up, we get 

J(t)>^{ [ (l-Mvj)ul(v)dy+ [ (l-fo,o(y))u 2 x (0,y)dy 

4 ^ J[T( y )>t]ns 3 v 7 i[r(j*)>t]n5 3 v 7 

+ / 4>iMul(y)dy+ / 2 ,o(y)«|(O,y)dy}, 

V[f(j/)<t<T(j/)]ns 3 ^[T(y)<«T(j;)]nS3 J 

since the last term on the right is zero as 53 C [T(y) > t]. A similar relation holds with S4 instead 
of 53. Consequently, putting together all this information about the various inequalities valid on 
the disjoint sets Si, S2, S3, and S4, we conclude by (3.44) that 

^-u(t,Z(t,y ))- ±u(t,£(t,y )) <-J (t) + ^. (3.45) 
at at kqtt 

To obtain now a suitable estimate on 
hmsup J (t + h) -J(t) = limgup I E(t + h,y)-E(t,y) ^ <{y) ^ + ^ ^ ^ 
HO n hio Js(t) n 

we distinguish two cases. If E(t,y) is the second component E 2 (t,y) of the minimum in (3.30), 
then by (3.37) and (3.41) we can estimate the contribution of the first term in (3.46) by zero from 
above. On the other hand, if the minimum is E x {t, y), then the first integral term in (3.46) is not 
larger (pointwise) than the nonnegative expression 



( E (t,y) + -^--J (t))<l> 1 ,o(y)ul(y) 
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in view of the estimates (3.40), (3.42), and (3.45). We conclude that 

l imsup J (* ± Q ~ J ^ < j( t ) + _ j ( t )) ||tz2|| L + K(j J (t). 

Since Jo(i) < 0, choosing the constant k = ||^x|| L i^ we now have 

Optimizing over all triples (V>°, </>?, </>2) we conclude 

J(u(t), v(t)) < J(u(0), v(0)) e 2t , t > 0. (3.47) 

Summing up the considerations made above, we proved the following result. 

Theorem 2. The trajectories t u(t) of (1.1) constructed in Theorem 1 are locally Lipschitz 
continuous as maps from [0, oo) into the metric space X equipped with the distance functional J. 
Moreover, the distance between two trajectories is also locally Lipschitz continuous as a map from 
[0, oo) into X . 



4 - Concluding remarks 

The following example shows that, in some sense, our distance functional J in (3.11) is "sharp" . 
Indeed, the convergence of the initial data in L°°(iR) n L 1 (iR) together with the weak convergence 
of the derivatives u x and u%, in L 2 (iR) does not guarantee the convergence of the corresponding 
solutions at later times t > 0. 

Example 1. Consider the functions f,g : [0, 1] i-> [0, 1] defined as 

(l-2x if xG [0 1/21 (1-Sx if xe [0, 1/6], 

f(x) = I t 1/2 11 ^) = i V2 if x € [V6 , 1/2], 

1 " " fc ■ I 1 -t if xG [1/2, 1]. 



Observe that 



f f(x) dx = f g'(x) dx = -l, f [f'(x)] 2 dx= f [g'(x)] 2 dx = 2 . 

Jo Jo Jo Jo 



Next, consider the function 

h(x) 

and define the sequences of initial values 
u n (x) -- 



1 — \x\ if \x\ < 1, 

if \x\ > 1, 



h(x) if x(£ [0,1], 

h{i/n) + y{nx-i + l) if xG [i^,±] i = l,...,n, 
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v n (x) 



h(x) if a; £[0,1], 

h(i/n) + ±g(nx - i + 1) if s e , i] i = l,...,n. 



Letting rnoowe now have the strong convergence \\u n — Vn|| L00 (jj) — * 0- Moreover, by construc- 
tion it is easy to see that at each point x G [0, 1], 

lim / ( (u n ) x (y) - (v n ) x (y) ) = lim / ((u n ) 2 x (y) - (v n ) 2 x (y)) dy = 

so that in L 2 [0, 1] one has the weak convergence 

(«n)i - («n)i , (uj 2 . - (V n )l , 

since both sequences are bounded in L 2 [0, 1] and the previous observation identifies the zero func- 
tion as the only possible weak limit. However 

u(t) = lim u n (t) 7^ lim v n (t) = v(t) 

n — >oo n — >oo 

for every t G (2/3, 1), where T = 2/3 is the time at which the gradients of the functions v n blow 
up. The last assertion follows at once from (2.27). 



We also would like to highlight the importance of requiring that the transport map ip in 
(3.10) be monotone nondecreasing. If in (3.11) we were to take the minimization over all maps ip, 
not necessarily monotone, we would obtain the classical Kantorovich-Rubinstein distance between 
measures, which generates the weak topology on the space of bounded, positive measures [V]. By 
restricting ourselves to monotone nondecreasing maps tp, the corresponding distance functional 
generates a much stronger topology. 



Example 2. Consider the sequence of Lipschitz functions 

u m {x) = < 



{ 

x — (i — l)/m 
^ i/m — x 



if x i [0, 1] 

if (t- l)/m < x < (2i- l)/2m 
if (2t- l)/2m < x < i/m 



, m . 



In this = ±1 and arctanua; = ±7r/4. The corresponding measures fj, u defined at (3.9) 

converge weakly to the measure /j, on M 2 x [— vr/2 , tt/2] defined as 

fj,(A) = - measjx G [0, 1] ; (x, 0, 7r/4) G a| + - measjx G [0, 1] ; (x, 0, -7r/4) G a| . 

In particular, these measures form a Cauchy sequence in the Kantorovich metric. However, these 
same functions u m do not form a Cauchy sequence w.r.t. the distance J. Indeed, let m < n. 
Consider the open intervals 



1 2i 



m 



2m 



2i 



2m ' m 



where u™ 



takes the values +1 and —1, respectively. Define the intervals I™~ similarly. Now 
consider any transportation plan (i/>, <t>i, $2)1 with tp non-decreasing. For each i = 1, . . . ,m, call vi 
the number of distinct intervals I™ + which intersect the image ■0(I™' + ). Since tp is monotone, if 
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Ui > 2, this implies that the image tJj(I™ + ) entirely covers v ; — 1 distinct intervals I™ . Because 
u™ = 1 on P i n+ and it" = —1 on each F-~ , this accounts for a cost > {yi — l)/2n. Next, if 
v\ + . . . + v m = n* < n, there must be n — n* intervals J?,^ , . . . , I^ f „»_„\ which do not intersect 

any of the sets for i = 1, . . . , m. These intervals must be contained in the image of some 

Z™ -- , or in the image of the set ip(M\[0, 1]), where u m = 0. This accounts for a cost > (n — n*)/4n. 

The above argument shows that, for any m < n, the cost of any transportation plan is bounded 
below by 

i ( m m 1 i _ 

JC**.*)^, „») > i- • max 5> - 1) , n - ^ > i- • ^ . 

For any fixed m, the right hand side approaches 1/8 as n — > oo. Therefore, the above is not a 
Cauchy sequence, in our transportation metric. 
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